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Abstract We have shown previously (Bobylev et al 2011) that some of the 
stars in the Solar neighborhood today may have originated in the same star 
cluster as the Sun, and could thus be called Solar Siblings. In this work we 
investigate the sensitivity of this result to Galactic models and to parameters 
of these models, and also extend the sample of orbits. There are a number of 
good candidates for the Sibling category, but due to the long period of orbit 
evolution since the break-up of the birth cluster of the Sun, one can only attach 
probabilities of membership. We find that up to 10% (but more likely around 
1%) of the members of the Sun’s birth cluster could be still found within 100 
pc from the Sun today. 

Keywords Solar neighborhood • Stellar motion • Stellar clusters and 
associations • Spiral galaxies • Astrobiology 


1 Introduction 

Stars are generally born in clusters. Therefore it is likely that our Sun was also 
a member of a star cluster some 4 billion years ago after it had condensed from 
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an interstellar gas cloud (Adams 2010). By now the cluster has dispersed and 
its stars have spread out over a considerable region of the Galaxy (Portegies 
and Zwart 2009). 

There are several reasons to try to trace the other members of this cluster 
which we call Solar Siblings. Just to mention one, the other members of the 
cluster must have had planetary systems of their own. Rocks flung out from one 
planetary system would sometimes have landed on a planet of another system. 
This could have lead to exchange of material between the systems, and even 
to exchange of life between the planets of the different systems (Valtonen et 
al 2009). This may have important implications for the search of life in the 
Galaxy. If the origin of life is a very rare process, it is possible that life in 
our Galaxy is confined only to the planets of the Solar Siblings, outside our 
Solar System. When evidence of life in exoplanetary systems will be perhaps 
uncovered one day, it is important to determine if such systems belong to Solar 
Siblings, or if they are members of the general star field in our neighborhood. 
Life among Solar Siblings alone would speak for the rarity of the origin of 
life. In this paper we ask whether it is still possible to recognize the Solar 
Siblings after some twenty Galactic revolutions and plenty of unpredictable 
orbit evolution in the past 4 billion years. 

One way to search for Solar Siblings is to find stars of the same age and 
metallicity as the Sun. The Sun has an exceptional metallicity for our neigh¬ 
borhood (Wielen et al 1996), with only 25% of otherwise similar stars having 
metallicity as high or higher than the solar value. Thus metallicity studies of 
the nearby stars may provide a good set of Solar Sibling candidates. Batista 
and Fernandes (2012) emphasize the chemical similarities and age in the search 
of the solar neighborhood, and pick out three candidate stars. Bobylev et al 
(2011) look at the orbital convergence of stars in the past in a model Galaxy 
with spiral arms and identify two Solar Sibling candidates. The star HD83432 
is common to both lists. Brown et al (2010) find 6 candidates using a smooth 
Galaxy model without spiral arms. One of the stars in their list is common 
with Batista and Fernandes (2012). Ramirez et al. (2014) suggest another star, 
HD162826, to be a Solar Sibling from a study of 30 candidates. In this work 
we vary the spiral pattern of the Galaxy in order to see the sensitivity of the 
result to the parameters of the pattern. Rather than looking for particular 
Solar Sibling candidates, we ask how likely it is that there are still today true 
Siblings near to us. 

Previous studies suggest that we may find as many as 10% of the original 
Sun’s birth cluster members within 100 pc of us (Portegies and Zwart 2009), 
while others (Mishurov and Acharova 2011) put the number as low as below 
0.1%. It is an interesting question which number is closer to truth. There 
are still large uncertainties as to the exact model of our Galaxy (Antoja et al. 
2009, Gerhard 2011) as well as to the secular evolution of our Galaxy (Sellwood 
2013). Therefore we try to cover a representative range of Galactic models. 
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2 Construction of orbits 

We calculate the stellar orbits by solving the following system of equations 
of motion based on a realistic model of the Galactic gravitational potential 
(Fernandez et al 2008): 

( 1 ) 

i) = + 217oC) 

r--— 

where is the Galactic gravitational potential; the (5, -q, Q coordinate system 
with the center at the Sun rotates around the Galactic center with a constant 
angular velocity Cq, with the 77 , and C axes being directed toward the Galac¬ 
tic center, in the direction of Galactic rotation, and toward the Galactic North 
Pole, respectively; Rq is the Galactocentric distance of the Sun. 

We use the Fellhauer et al (2006) Galactic potential model. The system of 
equations (1) is solved numerically by the fourth-order Runge-Kutta method, 
with the time step of one million years. 

In the Fellhauer et al (2006) model, the Galactocentric distance of the Sun 
is taken to be Rq = 8.0 kpc and the circular velocity of the Sun around the 
Galactic center is Vq = |l7o|-Ro = 220 km s“^. The axisymmetric Galactic 
potential is represented as the sum of three components—the central bulge, 
the disk, and the halo: 


^ — *^halo T ^disk “t” ^bulge- (2) 

In this model, 

- the halo is represented by a potential dependent on the cylindrical Galac¬ 
tic coordinates R and Z as <l^haio{R, Z) = Vq\yi{ 1 + R^/(P + Z'^jd?), where 
Vq = 134 km s“^ and c? = 12 kpc; 

- the disk is represented by the potential from Miyamoto and Nagai (1975) 
as a function of the same coordinates: 'd>disk{R, Z) = —GMd{R^ + {b + (Z^ -|- 
c)^/^)^)“^/^, where the disk mass Md = 9.3 ■ 10^° Mq, b = 6.5 kpc, and 
c = 0.26 kpc; 

- the bulge is represented by the potential from Hernquist (1990): d>i,uige{R) = 
—GMb/{R + a), where the bulge Mb = 3.4 ■ 10^° Mq and a = 0.7 kpc. 

If the spiral density wave is taken into account (Lin and Shu 1964, Lin et 
al 1969), then the following term is added to the right-hand side of Eq. (2) 
(Fernandez et al 2008): 

<Psp{R,0,t) = Acos[m{^2pt - 9) + x(i?)], (3) 

where 

^ _ (Rol?o)^/rotanz 


m 
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x{R) = 


tani \Rq J 


+ Xq- 


Here, A is the amplitude of the spiral wave potential; fro is the ratio of the 
radial component of the perturbation from the spiral arms to the Galaxy’s 
total attraction; fip is the pattern speed of the wave; m is the number of 
spiral arms; i is the arm pitch angle, z < 0 for a winding pattern; % is the 
phase of the radial wave (the arm center then corresponds to x = 0°); and x© 
is the Sun’s phase in the spiral wave. 

The spiral wave parameters are not very well determined (a review of the 
problem can be found in Fernandez et al (2001) and Gerhard (2011)). The 
simplest model of a two-armed spiral pattern is commonly used, although, 
as analysis of the spatial distribution of young Galactic objects (young stars, 
star-forming regions, or hydrogen clouds) shows, both three- and four-armed 
patterns are possible (Russeil 2003, Englmaier et al 2008, Hon et al 2009, 
Bobylev and Bajkova 2014). More complex models are also known, for example, 
the kinematic model by Lepine et al (2001) that combines two- and four¬ 
armed spiral patterns rigidly rotating with an angular velocity close to 12o- 
Note also the spiral ring Galactic model of Melnik and Rautiainen (2009) that 
includes two outer rings elongated perpendicular and parallel to the bar, an 
inner ring elongated parallel to the bar, and two small fragments of spiral arms. 
As applied to the Galaxy, the theories of nonstationary spiral waves with a 
fairly short stationarity time (several 100 Myr), a variable rotation rate, and a 
variable number of arms are also considered (Sellwood and Binney 2002, Baba 
et al 2009). The currently available data do not yet allow any one of the listed 
models to be unequivocally chosen. Therefore, here we apply the models of 
both stationary and non-stationary spiral pattern with different numbers of 
spiral arms, but disregard the influence of the bar. 

We choose Qp from the range 8 — 35 km s“^ kpc“^ (Popova and Loktin 
2005, Naoz and Shaviv 2007, Gerhard 2011). The pitch angle z is known rela¬ 
tively well and is between —5° and —7° for the two-armed spiral pattern and 
between —10° and —14° for the four-armed spiral pattern. The amplitudes of 
the velocities of the perturbation from the spiral density wave are 5 — 10 km 
s“^ (Mishurov and Zenina 1999, Fernandez et al 2001, Bobylev and Bajkova 
2010 ) in both tangential and radial directions. 

The Sun’s phase in the wave x© has a rather large uncertainty. For example, 
according to Fernandez et al (2001), this angle lies within the range 284°—380°. 
Having analyzed the kinematics of open star clusters, Bobylev et al (2008) and 
Bobylev and Bajkova (2010) find x© close to —120°, while the data on masers 
yield an estimate of x© = —130° ± 10° (Bobylev and Bajkova 2010); here, the 
minus sign implies that the phase angle is measured from the center of the 
Carina-Sagittarius spiral arm. 

According to the classical approach in the linear density-wave theory (Yuan 
1969), the ratio fro lies within the range 0.04-0.07 and the most probable value 
is frO = 0.05. The upper limit fro = 0.07 is determined by the velocity dis¬ 
persion of young objects observed in the Galaxy (at fro = 0.07, the dispersion 
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must reach 25 km s which exceeds a typical observed value of 10-15 km 
s“^). We also study fro = 0.1, see Antoja et al. (2009) and references therein. 

We adopted the following optimum parameter values: for the two-armed 
spiral pattern (m = 2), the basic pattern speed of the spiral wave f2p = 20 
km s“^ kpc“^, the pitch angle i = —5°, the maximum ratio fro = 0.05, and 
the Sun’s phase in the wave y© = —120°. We also took into account the Sun’s 
displacement from the Galactic plane Zq = 17 pc (Joshi 2007) and used the 
present-day peculiar velocity of the Sun relative to the local standard of rest 
{Uq,Vq,Wq)lsr = (10,11,7) km s“^ (Schonrich et al 2010, Bobylev and 
Bajkova 2010). 

In addition, we studied the models of Sellwood (2010) and Gerhard (2011) 
with the following parameters: 

1 - model of Sellwood (m = 4, = 18 km s“^ kpc“^), 2 - model of 

Gerhard {m = 4, f2p = 25 km s“^ kpc“^), 3 - model of Sellwood (m = 2, 
f2p = 8.1 km s“^ kpc“^). 

In the case of four-armed model (m = 4) the pitch angle is taken as i = 
— 13°. Motivations for this value of pitch angle can be found in Bobylev and 
Bajkova (2014). From this and a kinematic analysis made earlier we find that 
the most probable number of spiral arms is m = 4. In the case of two-armed 
Sellwood’s model we took i = —6°. 

Our initial condition is such that one thousand stars are placed uniformly 
inside a disk of radius 10 pc with a velocity dispersion of 1 km/s in the two 
coordinates of the Galactic plane. They represent the birth cluster of the Sun 
in this two-dimensional simulation. Then their orbits are calculated forward for 
4 billion years in a given Galactic potential, with the intermediate result noted 
at 2 billion years of integration. The pattern speed of the wave is varied from 
one calculation to another: its values are 8.1, 10, 15, 18, 20 or 25 km/s/kpc. 
Also we study the case of no spiral pattern, as well as the birth cluster being 
in exact corotation with the spiral pattern [Qp = 200/8 km/s/kpc), and 5 
km/s/kpc off the corotation. The relative force of the spiral component is 5 or 
10%, the pitch angle is —5° or —6° for two-armed cases and —13° for 4-armed 
cases. The phase of the Sun in the spiral wave is —120°. The cases include the 
models of Sellwood (2010) (see also Hahn et al (2011)) with m = 2 and m = 4 
spiral arms, as well as Gerhard (2011) favoured model with m = 4. We take 
the spiral arm force component to be either constant or variable as shown in 
Figure 1. 

In the real Galaxy the density waves probably do not last more than about 
six revolutions, as they form and decay and reform over the twenty revolution 
periods (Sellwood and Binney 2002, Sellwood 2012, Sellwood 2013). For this 
reason we construct models where the spiral arm strength varies in this man¬ 
ner. In the real Galaxy the evolution of the Galactic potential is likely to be 
more complex (Sellwood 2013); we use these models to illustrate the sensitivity 
of the nearby Solar Sibling numbers with respect to these parameters. 
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Fig. 1 The relative strength of the spiral pattern as a function of time in the model of 
variable spiral arms. 


3 Results 

After the integration we look at the distribution of the Solar Sibling stars in 
the Galaxy (Figures 2(a) - 2(d)). Figures 2(a) and 2(b) compare the streams of 
stars after two different time spans, 2 Gyr and 4 Gyr, respectively. Results are 
shown both for the case of no spiral arms, and for constant spiral arms of two 
different pattern speeds. In Figure 2(c) we compare the effects of spiral arms 
which are in corotation with the initial cluster, or are off from the corotation 
by 5 km/s/kpc in either direction. In Figure 2(d) the spiral arm strength varies 
as illustrated in Figure 1. 

To take a closer look at the differences between the different cases, we 
identify the center of each stream and study the R — 100 pc, i? = 50 pc and 
i? = 25 pc surroundings around this center. This is not necessarily where the 
Sun would be in our model. However, it has been argued that the Sun was 
initially near the center of its birth cluster, and that it has not suffered close 
encounters with other stars (Adams 2010); thus it may still be not too far from 
the center of the stream. With this assumption, the number of stars within 
the specified circles represent the expected number of Solar Siblings today. 

In Tables 1-4 we give the mean value of the number of the Solar Siblings 
for different models. The central values for i, Qp and fro are listed in Golumn 
(1) of Tables 1-4 while for the central value for x© = —120°, if not otherwise 
stated. The parameters of the spiral wave i, flp^ fro, X& were either fixed (Ta¬ 
ble 3) or varied over a wide range in accordance with a Gaussian distribution, 
using the mean values equal to the nominal ones and the standard deviation cr 
= nominal value / 2. The values beyond Scr were excluded. In Table 4 only the 




Probabilities for Solar Siblings 


7 


initial value of fro was varied, around the central value of fro = 0.1. In addi¬ 
tion to a constant fro, the time-varying f^o was studied in all cases, using the 
stated maximum value. We have carried out the Monte-Carlo statistical sim¬ 
ulations using 100 random realizations for three models (Sellwood-1, Gerhard 
and Sellwood-2). 

The Solar Sibling mean values at R=100pc and after 4 billion years are 
typically 16 ± 6. We conclude that even though the spiral wave parameters 
may be uncertain to some degree, we expect Solar Siblings in all three models. 

The maximum number of Solar Siblings (100 out of 1000 stars in the clus¬ 
ter) is found in the case of Hp = 10 km/s/kpc when R = 100 pc at the time of 
4 Gyr. This is rather a special case since the Galactic potential has refocussed 
orbits at 4 Gyr which were already more spread out at 2 Gyr. Even though 
it is an unlikely situation to be realized in the Galaxy, this possibility is also 
good to keep in mind. More typically, there are about ten Solar Siblings within 
the R = 100 pc radius, and possibly none within the i? = 25 pc radius from 
the stream center. 

Figures 3(a) and 3(b) plot the distance of a representative star from the 
cluster center in the constant spiral arm model and in the variable arm model, 
respectively, in the basic fro = 0.05 model. The star makes repeated close 
approaches to the center. We may assume that the Sun is also in the central 
region of the cluster so that this represents also an approximation to the dis¬ 
tance of the star from the Sun. We see that stars which make close approaches 
to the Sun do exist, but there is a certain amount of sensitivity with regard to 
the exact model. This is seen from the two lines which represent our standard 
variable spiral arm model, and another model where the spiral arm strength 
has been increased by 2%. Figure 3(c) shows a similar difference where the 
same initial orbit in the models of Sellwood (2010) with m = 4 and m = 2 
is compared with Gerhard (2011) model of m = 4, using the variable arm 
strength in each case. The diagram 3(d) for the constant arm strength shows 
qualitatively similar differences between the three models. 


4 Conclusions 

We have found that the Portegies and Zwart (2009) estimate of 10 to 80 Solar 
Siblings within the 100 pc radius, for the same number of initial stars as in our 
calculation, is surprisingly good, even though this study neglected the effect of 
spiral arms. On the other hand, on the basis of Mishurov and Acharova (2011) 
we would expect mostly zeros in our tables, which is clearly not the case. Part 
of the explanation of the difference is likely to be in the amplitude of the 
spiral wave which is 10% in the Mishurov and Acharova (2011) model while 
we use both 5% and 10% in this work. However, even our 10% models do give 
Solar Siblings, typically of the order of 10 (Table 2). Mishurov and Acharova 
(2011) also have randomly picked the ^2p values in the range of 14.2 and 51.8 
km/s/kpc which is another difference in the Galaxy model as compared with 
ours and it may also contribute to the different result. 
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Table 1 The mean number of stars within distance R (in pc) from the cluster centre 
after 2 and 4 G yr in different models of constant or variable density wave obtained using 
100 random realizations in a Monte-Carlo simulation where the spiral wave parameters are 
distributed over a wide range (cjj = nomma^j/2, frO,nominal = 0.05). 


R=100 R=50 R=25 R=100 R=50 

[pc] [pc] [pc] [pc] [pc] 

2 G yr 4 G yr 


R=25 

[pc] 


Sellwood: 

m=4 

Qp = 18 ± 9 km/s/kpc 

i=-13±6.5° 32 13 5 16 6 2 

/^o = 0.05 ± 0.025 
X0 = -120 ±60° 

Gerhard; 

m=4 

= 25 ± 12.5 km/s/kpc 

i=-13±6.5° 44 16 6 16 6 2 

fro = 0.05 ± 0.025 
X0 = -120 ±60° 

Sellwood: 

m=2 

17p = 8.1 ± 4.05 km/s/kpc 

i=-6±3° 53 22 8 21 8 3 

fro = 0.05 ± 0.025 
X0 = -120 ±60° 

Sellwood: 

m=4 

Qp = 18 ± 9 km/s/kpc 

i=-13±6.5° 44 19 6 18 7 2 

/r 0 =var(max 0.05 ± 0.025) 

X0 = -120 ±60° 

Gerhard: 

m=4 

17p = 25 ± 12.5 km/s/kpc 

i=-13±6.5° 38 17 6 18 8 3 

/r 0 =var(max 0.05 ± 0.025) 

X0 = -120 ±60° 

Sellwood: 

m=2 

17p = 8.1 ± 4.05 km/s/kpc 

i = -6 ± 3° 55 24 9 33 15 6 

/r 0 =var(max 0.05 ± 0.025) 

X0 = -120 ±60° 


Therefore the search for Solar Siblings is a worth while task; there should 
be a reasonable number of them in our near Galactic neighborhood. However, 
our study also shows the sensitivity to Galactic models. It is advisable that the 
candidates for Solar Siblings should be studied in a range of Galactic Models 
to see how robust the orbital convergence between the Sun and the Sibling 
candidate is. 
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Table 2 The mean number of stars within distance R (in pc) from the cluster centre 
after 2 and 4 G yr in different models of constant or variable density wave obtained using 
100 random realizations in a Monte-Carlo simulation where the spiral wave parameters are 
distributed over a wide range {ai = nominalil2^ frO, nominal = 0.1). 


R =100 

[pc] 

R=50 R=25 R: 

[pc] [pc] 

2 G yr 

=100 R=50 

[pc] ]pc] 

4 G yr 

R=25 

[pc] 

Sellwood: 

m=4 

= 18 zb 9 km/s/kpc 
i= -13 zb 6.5° 
frO = 0.1 zb 0.05 

X 0 = -120 zb 60° 

Gerhard: 

18 

8 3 

7 3 

1 

m=4 

Hp = 25 zb 12.5 km/s/kpc 
i = -13 zb 6.5° 
frO = 0.1 zb 0.05 

X 0 = -120 zb 60° 

Sellwood: 

25 

10 4 

10 4 

1 

m =2 

Hp = 8.1 zb 4.05 km/s/kpc 
i = -6zb3° 
frO = 0.1 zb 0.05 

X 0 = -120 zb 60° 

Sellwood: 

23 

9 3 

14 5 

2 

m=4 

Hp = 18 zb 9 km/s/kpc 
i = -13 zb 6.5° 
/■y.Q=var(max 0.1 zb 0.05) 

X 0 = -120 zb 60° 

Gerhard: 

19 

8 3 

8 3 

1 

m=4 

Qp = 25 zb 12.5 km/s/kpc 
i = -13 zb 6.5° 
/■^Q=var(max 0.1 zb 0.05) 

X 0 = -120 zb 60° 

Sellwood: 

27 

11 4 

10 4 

1 

m =2 

Qp = 8.1 zb 4.05 km/s/kpc 

2 = -6 zb 3° 

/■^Q=var(max 0.1 zb 0.05) 

X 0 = -120 zb 60° 

31 

12 4 

10 4 

1 
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Table 3 The mean number of stars within distance R (in pc) from the cluster centre after 
2 and 4 G yr in different models of constant or variable density wave frO = 0.05, 0.1. 



R=100 

[pc] 

R=50 

[pc] 

2 G yr 

R=25 

[pc] 

R=100 

[pc] 

R=50 

[pc] 

4 G yr 

II 

Sellwood: 

m=4 

Gp = 18 km/s/kpc 
frO = 0.05 

Gerhard: 

105 

61 

27 

81 

31 

8 

m=4 

Gp = 25 km/s/kpc 
frO = 0.05 

Sellwood: 

30 

10 

5 

16 

8 

4 

m=2 

Qp = 8.1 km/s/kpc 
frO = 0.05 

Sellwood: 

135 

75 

41 

66 

24 

7 

m=4 

Gp = 18 km/s/kpc 
y‘^Q=var(max 0.05) 
Gerhard: 

56 

17 

6 

25 

8 

1 

m=4 

Gp = 25 km/s/kpc 
/‘^Q=var(max 0.05) 
Sellwood: 

22 

4 

1 

13 

8 

3 

m=2 

Gp = 8.1 km/s/kpc 
/‘^Q=var(max 0.05) 
Sellwood: 

124 

74 

33 

62 

32 

19 

m=4 

Gp = 18 km/s/kpc 

frO = 0.1 

Gerhard: 

0 

0 

0 

0 

0 

0 

m=4 

Gp = 25 km/s/kpc 

frO = 0.1 

Sellwood: 

11 

6 

1 

13 

2 

1 

m=2 

Gp = 8.1 km/s/kpc 

frO = 0.1 

Sellwood: 

6 

0 

0 

3 

0 

0 

m=4 

Gp = 18 km/s/kpc 
/■^ 0 =var(max 0.1) 
Gerhard: 

17 

7 

0 

0 

0 

0 

m=4 

Gp = 25 km/s/kpc 
/■^ 0 =var(max 0.1) 
Sellwood: 

1 

0 

0 

2 

1 

1 

m=2 

Gp = 8.1 km/s/kpc 
y‘^Q=var(max 0.1) 

18 

8 

2 

16 

9 

4 
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Table 4 The mean number of stars within distance R (in pc) from the cluster centre after 
2 and 4 G yr in different models of constant or variable density wave obtained using 100 
random realizations in a Monte-Carlo simulation where /^o=0.1 di 0.05, other spiral wave 
parameters are fixed. 



R=100 

[pc] 

R=50 

[pc] 

2 G yr 

R=25 

[pc] 

R=100 

[pc] 

R=50 

[pc] 

4 G yr 

R=25 

[pc] 

Sellwood: 

m=4 

Qp = 18 km/s/kpc 
i = -13° 

21 

9 

3 

9 

4 

1 

frO = 1.0 zb 0.05 

Xo = -120° 

Gerhard: 

m=4 

Qp = 25 km/s/kpc 
i = -13° 

22 

9 

4 

14 

6 

2 

frO = 1.0 zb 0.05 

XO = -120° 

Sellwood: 

m=2 

Gp = 8.1 km/s/kpc 
i = -6° 

21 

9 

3 

9 

3 

1 

frO = 1.0 zb 0.05 

XO = -120° 

Sellwood: 

m=4 

Qp = 18 km/s/kpc 
i = -13° 

51 

22 

8 

7 

2 

1 

/■^Q=var(max 0.1 zb 0.05) 
XO = -120° 

Gerhard: 

m=4 

Qp = 25 km/s/kpc 
i = -13° 

12 

5 

2 

4 

2 

1 

/■^Q=var(max 0.1 zb 0.05) 
XO = -120° 

Sellwood: 

m=2 

Gp = 8.1 km/s/kpc 
i= -6° 

26 

11 

4 

17 

7 

2 

/■^Q=var(max 0.1 zb 0.05) 
XO = -120° 
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(b) 



(d) 


Fig. 2 The 1000 stars originate from a cluster marked by a large dot at time £ = 0. After 
2 Gyr panel (a) and after 4 Gyr (panels (b) - (d)) the stars have spread out in streams 
shown by small dots. Several cases of spiral pattern speed are shown as indicated in the 
figures. The center of each stream is also indicated. The circular orbit of the cluster center 
is depicted by a dotted line. In panel 2(d) the the spiral force component varies as shown in 
Figure 1, in other panels it is constant. 
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Fig. 3 The separation of a star from the cluster center assuming the constant spiral strength 
model (a) or variable spiral strength (b). The two lines show the result for two slightly dif¬ 
ferent values of the spiral potential strength. In panels (c) and (d) line 1 refers to Sellwood’s 
(2010) m = 4 model, line 2 to Gerhard’s (2011) model with m = 4 while line 3 refers to 
Sellwood’s model with m = 2. In panel (c) we use the variable spiral arm model, in panel 
(d) the constant arm model. 













